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We consider two species of particles performing random walks in 
a domain in R'* with reflecting boundary conditions, which annihilate 
on contact. In addition, there is a conservation law so that the total 
number of particles of each type is preserved: When the two particles 
of different species annihilate each other, particles of each species, 
chosen at random, give birth. We assume initially equal numbers of 
each species and show that the system has a diffusive scaling limit 
in which the densities of the two species are well approximated by 
the positive and negative parts of the solution of the heat equation 
normalized to have constant P norm. In particular, the higher Neu¬ 
mann eigenfunctions appear as asymptotically stable states at the 
diffusive time scale. 

1. Introduction. A branching particle system representation for the heat 
equation solution with positive temperature was introduced in [4] and later 
studied in [5] (see also [15] and [16]). Here is an informal description of 
that model and one of the main results, proved in [5]. Suppose that D is 
an open set in and N Brownian particles move independently inside D. 
Whenever one of these particles hits the boundary of D, it is killed and one 
of the other particles, randomly chosen, splits into two particles, so that the 
number N of particles remains constant. When the number of particles goes 
to infinity and the initial (normalized) distribution of particles converges to 
a measure on D, then the particle density converges to the solution of the 
heat equation for every time t >0, with the appropriate initial condition, 
normalized so that it has a constant total mass for all times t>0. 

The main purpose of this article is to study a related model that involves 
two different types of particles. We call them + and — particles. The two 
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types of particles perform independent symmetric random walks with re¬ 
flection on the boundary, and annihilate each other on contact. When two 
particles of different signs annihilate each other, then two other particles are 
chosen randomly, one -|- particle and one — particle, and each one of these 
particles splits into two particles of the same type as the parent particle, so 
that the numbers of -|- and — particles are conserved. 

The reader might have noticed that we have changed the Brownian par¬ 
ticles to symmetric random walks. This is because, except in dimension 1, 
the Brownian particles would not meet, so one would have to have them 
annihilate when they came within some e > 0 of each other. So we might as 
well just put the particles on a discrete lattice. Hence, particles in our new 
model will be represented by random walks on the lattice n D, with £~'^ 
comparable to the number of particles. We will assume that particles reflect 
from the boundary of D in a way to be described precisely later. 

One can easily generalize the model to have more than two types of 
particles. Particles of different types annihilate on collision and two other 
particles, from those two types, are chosen at random to split into new 
particles. The main effect in these systems is segregation of types. At a 
macroscopic scale one sees the region D decomposed into disjoint regions 
Di{t),... each occupied by one of the types. The density of type i 

evolves on Di by the heat equation with creation, with Dirichlet boundary 
conditions on boundaries internal to D, and Neumann boundary conditions 
on boundaries coincident with the boundary of D. The regions Di them¬ 
selves evolve according to the behavior of the type densities on either side 
of the boundary. In general, the boundary dynamics are not simple. At the 
present time the probabilistic picture is based purely on simulations and 
informal calculations. See [4, 9, 10] and [11] for details. One expects deter¬ 
ministic evolution on the macroscale, as described above. Note that related 
limiting equations for segregation can be derived in the scaling limit of ap¬ 
propriate systems of reaction-diffusion equations (see [6, 7] and references 
therein). These problems are quite difficult because, unlike most problems 
in hydrodynamic limits, the invariant measures of these systems are either 
unknown or not well understood. Of course, well inside the regions where 
each individual species resides, one expects a product of Poisson distribu¬ 
tions. However, it is not a priori clear that such regions exist, and anyway, 
all the interesting behavior in the system is at the boundaries between the 
species. 

In this article we consider the special case of two types, and assume, in 
addition, that the number of each type is the same, say, N, where N is 
proportional to the number of lattice sites, so that the densities are of or¬ 
der 1. Our main result is that, as N ^ oo, the densities of the two types 
of particles converge to the positive and negative parts of the solution p 
of the heat equation normalized to have the correct total masses. Letting 
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D^{t) and D^{t) denote the supports of positive and negative parts of the 
solution, we can see that this corresponds to the general picture just de¬ 
scribed, with a particular, nontrivial, evolution of the boundary between 
and D-. The evolution is in the normal direction u to the bound¬ 
ary and proportional to (A -|- V)p+/dpp+ = (A -|- V)p-/d^p- there. Here 
V = j£)S{p = 0)\Vp\‘^dx/ fjj \p\dx is the creation term at time t needed to 
keep the masses conserved. 

The method we use is not standard in hydrodynamic limits. The main 
methods known in hydrodynamic limits are based on entropy or attractive¬ 
ness. The model is not attractive, and in problems like this one, where the 
object of interest is a boundary separating two regions on which live mu¬ 
tually singular local equilibrium measures, entropy methods do not seem to 
be useful. So it is important to develop new methods that can be used for 
such problems. 

The proof is based on an observation that the microscopic evolution of the 
density rj of the difference between the occupation numbers of the two types 
of particles takes a particularly simple form dr] = {/S. + V)r]dt + dM, where 
M are martingales, A is the discrete Laplace operator, and V = V(t) is the 
rate of annihilation. Some calculus and the particular form of the equation 
allow us to control V. Because the proof depends on the particular form 
of the equation, it does not generalize in an obvious way to systems with 
more than two types or with unequal numbers of particles. However, it does 
suggest the precise form the result should take in other, more complicated, 
cases. 

2. Preliminaries. We collect in this section a few results that may be 
known but we could not find a ready reference for them. 

We start with a discussion of a single particle model—reflected random 
walk on a lattice region approximating a region D C M'^. Let D he a bounded 
connected open set in d>2, and let n D for e > 0. We want 

to construct a nearest-neighbor continuous time random walk Wf on D^. 
Moreover, we want VTf to converge to reflecting Brownian motion in D as 
e ^ 0. A natural choice for Wf would be a process that jumps to one of 
its neighbors with equal probabilities, after a suitable time delay. In the 
main part of our project, such a process is somewhat inconvenient to deal 
with. We will use a different model for a “reflecting random walk” on D^, 
described below. 

The process is a hnite state continuous time Markov process so it is 
fully described by the jump distribution and holding time at each state. We 
start by describing its jump distribution. Let dD^ be the set of x £ which 
have fewer than 2d neighbors in D^. If x £ D^\ dD^, then the random walk 
Wl will jump from x to any of its neighbors with equal probabilities. 
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We will assume from now on that D has an analytic boundary. This is used 
later to get some easy estimates on the lattice Laplacian of eigenfunctions 
of the true Laplacian. One expects the results of this article to hold with 
much weaker assumptions on the boundary, but we have not pursued this 
here. 

If e is very small and x G dD^ has fewer than d neighbors in D,., then this 
implies that the normal vector to the boundary close to x is almost parallel 
to one of the axes. It is not hard to see that by removing from all x G dD^ 
with fewer than d neighbors in D^, we obtain a set D'^ with the property 
that all X G dD'^ have at least d neighbors in D'^. By abuse of notation, we 
will refer to this set as D^. Note that every point in dD^ lies at a distance 
from dD not exceeding 2e. 

For X G dD^, let x' G dD be the closest point to x on dD and let n(x') be 

the inward unit normal vector at x'. Let {ei(x') = n(x'), 62 ( 3 ;'),... , 0 ^( 3 :')} 
be an orthonormal basis, depending on x. Let p^y denote the probability 
that makes the next jump from x to y (assuming it is at x now) and 
note that, by assumption, 

(2.1) E = 

y£Ds,\x-y\=e 

We want to find pxy for x G dD^ so that, for some ci > 0, 

(2.2) E Pxy{y-x) = cin{x'). 
yeDe,\x-y\=e 

Note that equations (2.2) impose only d — 1 constraints because ci > 0 is 
arbitrary. Hence, (2.1) and (2.2) effectively form a set of only d equations, 
and we have at least d unknowns pxy It is not hard to see that (2.1)-(2.2) 
have a solution. 

Next we choose the holding times. For x G De \ dD^, we let the hold¬ 
ing time at x have the mean he{x) = . This corresponds to the usual 

space-time scaling for Brownian motion. For x G dD^, the mean h^{x) of 
the holding time at x is chosen so that 

(2.3) hmi E E{{(Wt^,-Wt)-By{x')f\Wt = x) = l. 

^ l<3<d 

The discrete Laplacian A^, that is, the generator of the process is given 
by 

(2.4) Aef{x) = h-\x) E Pxyifiv)- f{x)). 

y&De,\y-x\=e 

The following lemma provides an intuitive basis for the interpretation of 
our main results. Since the lemma is not used later in the paper, we only 
observe that it can be derived from Theorem 6.3 of [21] and we do not supply 
a formal proof. 
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Lemma 2.1. Suppose that G and ^ xq & D as e — > 0. If Wq = 
Xs, then {Wf,t>0} converge weakly to the reflecting Brownian motion on 
D starting from xq, as s ^0. 

We have assumed that D has an analytic boundary. This implies that D 
has a discrete spectrum for the Laplacian with Neumann boundary condi¬ 
tions; see, for example, Section 2 of [1]. In the Dirichlet case, the spectrum 
is discrete for any bounded open connected set D. In the Neumann case, the 
spectrum is not necessarily discrete if D is bounded but has a nonsmooth 
boundary; see, for example, [17]. It is quite easy to see that the spectrum 
is discrete for the Neumann Laplacian when dD can be represented locally 
by a Lipschitz function; see, for example. Section 2 of [1]. This condition 
on the boundary can be substantially relaxed—it is enough to suppose that 
the boundary is locally represented by a function which is Holder contin¬ 
uous with sufficiently small exponent, see, for example, [19]. As we said, 
our assumptions on the smoothness of dD are much stronger than that, for 
technical reasons. 

Let —An, re = 0,1,2,..., be the eigenvalues of the Laplacian on D with the 
Neumann boundary conditions, and let fn be the corresponding eigenfunc¬ 
tions. We list the same eigenvalue more than once, if it has a multiplicity 
greater than 1. For a measure /i on D, we let jin = J-jy(t>ndp,. 

The following review of random measures is based on Section 3 of [12]. 
A sequence of finite measures fD on a space A converges weakly to p, if for 
any bounded continuous function /, we have f{x) dgl'jx) fj^ f{x) dfj,{x) 
as re —> oo. Let MpjA) denote the space of all hnite measures on A equipped 
with the topology of weak convergence, let Mi (A) be the subspace of MpjA) 
consisting of probability measures, and let Mp^flA) denote the subspace of 
Alp (A) consisting of measures with total variation less than or equal to c. 
Let B{D) denote the family of Borel subsets of D. Suppose that, for some 
probability space (12, D, P) , the function ji: 13(D) x 12 ^ M has the following 
two properties: (i) for a fixed to, is a finite measure on D, and (ii) 

for a fixed A € 13(D), iJ,(A,-) is a random variable. Then the distribution 
of pL is an element of Mi(Mp(D)). We will refer to ^ as a “random mea¬ 
sure.” If A is a metric space, then the space of right-continuous functions 
with left limits, /: (0, oo) —> A, equipped with the Skorohod topology (see 
Section 3.6 of [12]) will be denoted by 5((0,oo),A). We will use the space 
AIi(S((0,oo), Mp(D))) to state our main results. To give a meaning to this 
symbol, we have to specify a metric on Mp(D). We will use the Prohorov 
metric which induces a topology on Mp(D) equivalent to weak convergence 
of measures (see [13], Chapter 3, Sections 1 and 3). In this article we will 
be concerned with the convergence of processes on the open half-line (0, oo), 
not the usual semi-closed half-line [0, oo). The convergence in the Skorohod 
topology on (0, oo) is defined as the convergence in the Skorohod topology 
on every compact subinterval of (0,oo). 
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Lemma 2.2. Suppose that f :D is a continuous (and, hence, bounded) 
function and let an = Jjy f ix)(j)n{x) dx. Suppose that ci < oo and pL^, k>l, 
are random measures (possibly defined on different probability spaces), with 
distributions in Mi{Mp^ci{D)) ■ Assume that, for every fixed n, limk^oo fin — 
On in distribution. Then the distributions of pf converge weakly in M\{Mp{D)) 
to 5fj_, where fj,{dx) = f{x)dx. 

Proof. Let dv^{x) = dfj.^{x) — f{x)dx and note that ^ 0 in distri¬ 
bution as k ^ oo, for every n. It will suffice to show that the distributions 
of converge weakly in Mi{Mp{D)) to do, where 0 is the measure identi¬ 
cally equal to 0. Since D is compact, so is Mp^cfiD) and it follows that the 
sequence is tight and contains a convergent subsequence. Let n be the 
weak limit of a subsequence of . It will be enough to show that = do- 

Let Gn'.Mp{D) —> M be defined by Gn{cr) = \ Jjycpnix) da{x)\ A I. The 
functionals Gn are continuous and bounded, and distributions of con¬ 
verge weakly to u in Mi{Mp{D)), so for every fixed n, EGn{n^^) —> EGn{n) 
as j —> oo. By assumption, z>^ ^ 0 in distribution for every n, so EGni^) — 
EGn{n^^) = 0 for every n. Hence, n is supported on measures a G 
Mp{D) with the property that an = 0 for all n. It will suffice to show that 
every measure with this property and finite total variation is identically 
equal to 0. 

Fix a nonrandom measure a on D with a finite total variation and such 
that cJn = 0 for all n. Fix any Borel set Ad D. li will be enough to show 
that JjylA{x) da{x) = 0. 

According to the Weyl formula, An ~ (see [8], Vol. I, Chapter VI, 
Section 4.4; or see [19] for a recent strong version of this theorem). By 
Theorem 1 of [14], ||())n||oo < so ||(/>n||oo < _ 

Let Pt be the transition semigroup for the reflected Brownian motion in D 
and for some fixed t > 0, let g{x) = Pt'^A^x). Then g(x) = Y)n{^A)n4>n{x)e~^‘^^, 
where (iA)n = f-^lA{x)(t)n{x) dx. If we write gix) = Jfn dnfinix), then 

\9n\ < \D\Un\\oce-^-^ < 

Without loss of generality, assume that the total variation of a is not 
greater than 1 and note that 

[_\4>n{x)\d\a{x)\ < |D|||(/)n||oo < 

JD 

This and other estimates obtained so far imply that 

n. ^ n 
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Hence, we can change the order of integration and summation in the follow¬ 
ing formula: 

_g{x)da{x)= /_ dfT(x) = _(j)n{x)da{x) = 0. 

Jd Jd ^ „ Jd 

We have proved that fjyFtlA(x) da(x) = 0 for every t > 0. Clearly, \PtlA{x) \ < 
1 for all X, and PtlAix) —> 1a( a^) for almost every x G H, so by the domi¬ 
nated convergence. 


fiAix) da{x) = lim / da{x) = 0. 

Jd 


□ 


The following result can be proved using standard methods, so we leave 
its proof to the reader. 

Lemma 2.3. Suppose that, for every e > 0, we have a real-valued process 
{R^{t),t > 0} which is equal to i?f (t) -|- + Rsit), and these processes 

satisfy the following conditions: 

(i) For every e > 0, t^ Rfit) is right-continuous and nondecreasing a.s. 

(ii) For every fixed t>0, the family {i?f (t),e > 0} is tight. 

(hi) For every fixed t > 0, supo<s<t |1^2('S)I converges to 0 in distribution, 
as e —i- 0. 

(iv) For every t>0. 


\RUt2)-RUti)\>S2'^=0- 

Then the family {R’^,£ > 0} is tight in Mi(S'((0, oo),] 


lim sup lim sup P ( sup 

< 5 i , 52—>0 £—>0 \0<ti,t2<t,\ti—t2\<6i 


3. Main results. We will now describe the main object of our study, a 
particle system. Our description will be partly informal because this will 
make it more accessible without loss of rigor. 

The state of the particle system at time t > 0 will be encoded as an 
integer-valued random function g{t) = gx{t) on D^. We will often suppress t 
in the notation and write px, with the convention that if gx > 0, then there 
are gx particles of type + at x € D^, and gx <0 signifies the presence of \gx\ 
particles of type — at x. Obviously, gx = 0 means that there are no particles 
at X. Since -|- and — particles annihilate each other instantaneously, we do 
not have to have a notation representing both types of particles at the same 
site of D^. We assume that there are N particles of type -|- and N particles 
of type — for every t. 

The easiest way to describe the evolution of g is to use the particle picture. 
Each particle performs continuous time symmetric simple random walk (de¬ 
fined in Section 2), independent of other particles, until one of the particles 
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hits a particle of the other type. When a particle jumps to a site occupied 
by at least one particle of the different sign, two particles of opposite signs 
annihilate each other. At the same time, one particle is chosen uniformly 
from the family of + particles and one — particle is chosen uniformly as 
well. Each of the two chosen particles splits into two offspring of the same 
type as the parent, so that the number of particles of each type remains 
constant. The particles move independently until the next annihilation and 
birth event, and the evolution continues in the same manner. From the point 
of view of the mathematical description of the model, it is more convenient 
to represent the “annihilation and birth” event as a jump of the annihilated 
+ particle to a randomly chosen + particle, and the same for the annihilated 
— particle. 

In the above description of the dynamics of the system, when we say 
that a particle is chosen “uniformly,” it means that we choose one of the 
N particles with the same probability 1/A^; in other words, a particle that 
is annihilated may be the one that splits into two offspring. In such a case, 
the offspring are born at the site where the annihilated particle resided just 
before the jump that lead to its annihilation. This is a different convention 
than in [5], but this convention will simplify some formulas and, of course, 
it makes little difference when N is large. 

The informal description given above can be rigorously expressed in terms 
of the generator L = for the process. Let = max(a,0) and a~ = 
max(—a,0). The configuration rj such that r}x = l and r/y = 0 for y ^ x will 
be denoted \x- The formula for the generator L of the process y is 

Lfiv) 

= K^ix)Pxy\vx'^{vy>o}{f{Tl-'^x + ly)-f{r])) 

x,y£Ds I 

+ l{77^<0}(/(?/ + Ix - ly) - fiv)) 

+ 'nt^{Vy<0}^ ^ 

X ^^%(/(^-Ix + Iy +In-I d) -/(?/)) 

X VuVvifiV + '^x-'i-y-lu + 'i-v)-fiv)) 

We would like to point out that the first sum in the above formula extends 
over only those G that satisfy \x — y\= e. This extra condition is 
enforced automatically by the presence of the factor pxy. A similar remark 
applies to other formulas in this paper. 
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The normalized particle density u{x,t) = is defined by u{x,t) = 

N~^e~'^rix for x & D^. Typically, we will be interested in the population size 
of order N Let Pt denote the heat semigroup on D, corresponding to 
the reflecting Brownian motion, and for a measure /r on D, let PtH denote the 
measure with the density JjyPt{x, ■) dfi{x). When fJ- = 
and \x denotes the probability measure with the unit mass at x (by abuse 
of notation), then we will write Ptu^’’^{0) to denote PtiJ.- In other words, 
dy = Y.x&De 0)P^x, y) dy. 

For a nonzero measure /r on T) of finite total variation, we define JI to 
be cy, where c = c(/r) is a normalizing constant chosen so that the total 
variation of JI is equal to 2. For definiteness, we let /I = 0 if /r = 0. 

For the meaning of Mi(5((0, oo), Mf(Z)))), see Section 2. 

Theorem 3.1. Suppose that e —> 0 and N ^ oo in such a way that 
ci£~^ < N < C 2 £~'^ for some constants 0 < ci < C 2 < oo. Assume that D 
has an analytic boundary, u^’^{x,Q) are nonrandom, and for some signed 
measure y on D which does not vanish identically, J2x&De O)!^; —>■ y 

in Mf{D) as e —>■ 0. Then, as e —!• 0, 

^ e'^u^’'^{x,-)\x-Pp^do in Mi{S{{Q,oo),Mf{D))). 

x£De 


The symbol do in the last formula refers to the process identically equal 
to the null measure. 

Note that we have not assumed that the total variation of y is equal to 
2. The total variation of y cannot be larger than 2, but it can be smaller 
than 2 when the particles are tightly interspersed at time 0. If the total 
variation of y is less than 2, then for small e, the particle configuration has 
an almost instantaneous jump at time 0 to a configuration that approximates 
JI. For this reason, we obtain convergence only in Mi{S{{0,oo), Mf{D))), 
not in Mi{S{[0,oo), Mf{D))). If the total variation of y is equal to 2, our 
arguments show that the convergence holds in Mi(S'([0, oo), Mi7’(T)))). 

The assumption that u^’^{x, 0) are nonrandom measures is easy to remove— 
we added it for technical convenience only. 

Theorem 3.1 is a special case of Theorem 3.2 below. We need some more 
notation to present this more general result. 

Recall from Section 2 that —A^, n = 0,1,2,..., are the eigenvalues of 
the Laplacian on D with the Neumann boundary conditions, and (pn are 
the corresponding eigenfunctions. For a measure y on D, we write JJn = 
f-^pndy. For a function /:T»£ we let {f,g) = e'^Y.xeD, fix)g{x). The 
Fourier coefficients for the “density” u{x,t) will be denoted Un = Un{t) = 
{u{t),(j)n)- In other words, Un is the nth Fourier coefficient for the measure 
J2x£D, e‘^u^’%x,t)ix. 
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Note that if J2x&De —> At in Mp{D) as e —> 0, where /r is a 

signed measure that is not identically equal to 0, then, for some n, there 
exists a > 0 such that |h„(0)| > a for sufficiently small e > 0 (see Lemma 2.2 
and its proof). 

Theorem 3.2. Suppose that e —> 0 and N ^ oo in such a way that 
ci£~‘^ < N < C 2 £~'^ for some constants 0 < ci < C 2 < oo. Assume that D 
has an analytic boundary and for some n, there exists a > 0 such that 
infjv,e |h„(0)| = a. Then, as e —> 0, 

^ Oix - T.u^’^(O) ^ (5o in Mi{S{{0,oo),Mf(D))). 

x£De 

The remainder of the paper is devoted to the proof of this theorem. We 
start with a very informal overview. The proof will be divided into several 
steps. 

In Step 1 we observe that the density field, u, of + particles minus — 
particles, satishes a particularly simple, linear equation at the microscopic 
scale, 

du, = [A: + V,]u, + dM„ 

where A* is the adjoint of the random walk generator, V = V{t) is the 
particle annihilation rate, and M is a field of martingales. 

This is based on the following elementary observation. Let A be the gen¬ 
erator of a continuous time Markov process with state space S. Let L denote 
the generator of a system of particles of two types, on S, performing this dy¬ 
namics, and, in addition, annihilating on contact. Let r]x denote the number 
of the first type minus the number of the second type, at x G S'. Then 

L{r],i:) = {A*rj)a:. 

Adding the particle creation term, which is clearly of mean field type, gives 
the preceding linear equation. 

Now our problem is to show that ne(t,x) is close to v{t,x) = 
v{t,x) satishes a similar looking equation dtv = [A -|- W]v, where the job 
of IT = W{t) is to maintain the total mass of v. Formally differentiating 
fj^\v\dx gives W{t) = fpd(v = 0)|Vnpdx/ fj^\v\dx. So there are two key 
things to prove: First, that the martingale terms vanish in the limit. Second, 
that T looks like IT on the macroscale. Note that control of Mg itself is not 

enough. One really needs to control martingale terms like /g eT dMf:{s) 

The T are the main unknowns, and this is the key difficulty of the problem. 

The linear form of the equation suggests the use of Fourier analysis. In 
Step 2 we rewrite the equation as a system of stochastic differential equations 
for the Fourier coefficients of the density held [see (3.6) and (3.7)]. A simple 
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observation is that, because of the conservation law, and well-known bounds 
on the eigenfunctions of the Laplacian, the Fourier coefficients themselves 
are bounded. A bit of calculus, and the particular form of the system, then 
allow us to obtain an a priori estimate on 14 [see (3.15)]. Once this is done, 
a preliminary form of the limiting equation can be obtained [see (3.16)]. 

(3.16) tells us that what we see macroscopically at time t is the evolution 
of the density field by the heat equation, as long as we are willing to multiply 
by some scaling factor. At first glance, this would appear to imply the full 
result, because the scaling factor is fixed by the conservation law. However, 
there is still a lot of work to show that this is in fact the case. What could be 
happening is that the two type of particles are coexisting on some mesoscopic 
scale. All we would see on the macroscale is a net decrease in the total 
mass. In Step 3 we show that this cannot happen. The reason is that such a 
situation would lead to a very large rate of annihilation, as -|- and — particles 
would find themselves unnaturally close to each other. So V) would get very 
large. But we already have a bound (3.15) which prevents this. Finally, in 
Steps 4 and 5 similar ideas are used to show the required tightness. 

We now proceed with the proof. 


Proof of Theorem 3.2. Step 1. In this step we will show that, for 
zGD^ and /(ry) = 


(3.1) 


LJ{rj) = + Pry, 


where V = V{e, N, ry) is the (normalized, instantaneous) jump intensity for 
the particle system in state ry, defined by 

V = 2N-^ hj^ix)Pxyir]t^7jy<0} + VxMvy>0}) 

x,y£De 

and 

Kfix) = {h-^{y)pyxf{y) - hj\x)pxyf{x)). 

y&De 

This operator is the adjoint of the discrete Laplacian Ag given by 
A£/(x) = h-^{x) Y Pxyifiv) - fix))- 

y&De 

Fix some z G and let /(ry) = ry^. Then 

Lf{v)= Y ??^l{^^>o}(l{y=H - 1 h=4) 

x,v£Dp_ I 


~^Px'^{riy<0}i'^{x=z} l{j/=2}) 
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(3.2) 


The sum 


^ 7 y Vu Vv i^{y=z} T ^{u=z} 

u^vGDs 

^{x=z} ^{v=z}) 

+ 'nx^{rty>0}^ ^ 

^ ^ i^{x = z} T ^{v=z} 

u^vGDs 

— l{y=2} “ l{u=2}) j 

of the terms in (3.2) with the indicators 1{3;=2} is equal to 

-'nt^{r!y>0}+Vz^{r!ym 
-vt^{r,y<0}N~‘^ 

u,vGD£ 


y&De 


'zy 


(3.3) 


+ Vz^{rjy>0}N ^ J2 PuVv 

u,vGD£ 

= J2 K^i^)p^yi-vt'^{vym + Pz^{vym 

y&De 

- r}tHyy<0}N-^N^ + il-l{r,y>0}N-‘^N‘^] 

= - X! K^i^)PzyVz- 

y&De 

A similar calculation shows that the sum of the terms in (3.2) with the 
indicators l{y= 2 } is equal to 

(3.4) ^ h~^{x)pxzrix- 

x&De 

The sum 


x&De 

of the terms in (3.2) with the indicators '\-{u=z} is equal to 


{x)Pxy 

x,y£De 


Px'^{yy<0}N ^'ntPv 

v&De 

- Px^{ny>0}^~‘^ Y1 ^zVv 

v&De 

-- h-^{x)p:,y[ritl{r,y<o}N~^ptN-r}-l{r,y>Q}N~^V'^N] 

x,y&De 


(3.5) 
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= 

Similarly, the sum of the terms in (3.2) with the indicators 1(^=2} is equal 
to ( 1 / 2 ) 1 / 772 . Combining this and (3.2)-(3.5), we obtain (3.1). 

Ste'p 2. We will now derive estimates for the Fourier coefficients of u. 
Note that 77 is a finite state continuous time Markov process. Therefore 
(see [13], Chapter 4, Proposition 1.7 and Chapter 4, Section 2), M„(t) = 
Unit) — jQLun{s)ds is a martingale. In other words, the Fourier coefficient 
Un satishes the following stochastic differential equation; 

(3.6) dun = Lun dt + dMn- 

Since (/>„ is the nth. eigenfunction, = —Xn(f^n, where A is the Laplacian 
on D with the Neumann boundary conditions. Recall that we have assumed 
that D has an analytic boundary. It was pointed out in [2] that, in view 
of Theorem 5.7.1 on page 169 of [18], if the coefficients of a second-order 
elliptic equation are real analytic on a bounded analytic domain D up to 
the boundary, then, for every point 2 : on dD, there exists a ball B centered 
at 2 ; such that solutions of the elliptic equation can be extended to be real 
analytic functions on B. Hence, for every point x € D, (pnix) is equal to its 
power series in a neighborhood of x. The radius of convergence of the series 
is strictly positive for every x € D and continuous as a function of x. Since 
the domain D is bounded, the radius of convergence is bounded below for 
all X G H by a strictly positive constant. The coefficients of the power series 
are continuous as functions of x. We will estimate J2yeDe Pxyi4’niy) — 4’nix)) 
using power series expansion at x. Since the normal derivative of 4>n vanishes 
on the boundary, (2.2) implies that the sum of linear terms vanishes. The 
sum of quadratic terms is equal to e^A(/„(x) + o(e^) because of (2.3). The 
contribution of higher terms is o(e^). Hence, the sum in question is equal to 
s‘^A(j)n{x) +o(e^). It is easy to check using the dehnition (2.3) that h^{x) is 
of order e^. All these estimates imply that we have 

A^Cpnix) = A4>nix) + Ipeix) = -Xn<Pnix) + '<Peix), 
where |7/e(3^)] < cin)e. We have Lr]x = A*ijx + Vilx and 

Un = {u,4>n) =e'^ ^ N~^£~'^'qn,(t)nix) = N~^ ^ ?7a:</n(a^), 
x^De x£Ds 

so 

LUn = N~^ ^ iLr]^)4>nix) 

xSDe 

= N~'^ {Alr]^)(j)nix) + N~^ Y iVr]x)4'nix) 
xSDe x&De 
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= ^ ^ I] ( H ^iy)PyxVy - K ^ix)PxyVx) I 4>n{x) 
xeDg \yeDe / 

+ N~^V rjxMx) 

x&De 

(3.7) =N~^ (h~\x)pa;y((>n(y) - h~^(x)pa;y((>n(x))Ua; 

x&De \yeDe / 

+ 1V“V Y VxMx) 

x&De 

= 1V^ Y 71x(^e(pn{x)) + VN~^ Y ^xf^nix) 

xSDe x£De 

= N~^ Y rix{-K4>n{x) +ll^e(x)) + VN~^ Y 'Hx^’nix) 

x£De xGDs 

— '^n)'^n “t“ 

where |'he,n| < c(n)£. 

The process p has only a finite number of states, so Lun is uniformly 
bounded. The process Un jumps after an exponential waiting time. These 
facts and the formula Mn(t) = Un(t) — Jq Lun(s) ds imply that 

E[(Mn(t + s)-Mn(t))^] 

™0 E[iUn{t + s)-Unim 

It follows that 

(3.8) E[M^it)]= ['E[Anis)]ds, 

Jo 

where An = An(t,p) is given by the following formula: 

= lim(l/s).E[(u„(t + S) - Un(t))‘^\p{t)] 
s^O 

x,yeDe I 

+ l{»?^<O}(0n(a:) - </>n(y))^ 

+ vt^{Vy<0}]^~‘^ Y vtVviMy)-<Pn(x) 

u^vGDs 

+ 4'n(u) -(t)n{v)f 

+ dx'^{ny>0}N~‘^ Y yudt{(l>n{x)-(i)n{y) 

U^VGDs 
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- (f)n{u) + 4 )n{v)Y 


This implies the following bound for for small e: 

(3.9) An < CiN-^V{s^\\V(^n\L + UnWL) < PnN-^V, 

where jSn < oo depends on (j)n, and V stands for the usual gradient acting 
on functions defined on It follows from (3.6) and (3.7) that 

n„(t) = +^%-/o 

(3.10) =h„(0)e/o(^M-^")'^'-+ 

Jo 

Jo 

= h„(0)e/o(^M-^'^)'''- + Rn^Yt) + Rn,2{t). 

n ^ 

Since Rn,i{t) = /q eT dM^, (3.8) and [20], Chapter II, Section 6, 

Corollary 3 and Theorem 29, show that 


E[RlYt)]=E 




In view of (3.9), 


E[RYYt)]<(JnN-^E 


uo 




We have 


Jo 

= \n + (l/2)(e2/o('^W-^")'^" - 1) 

Jo 

^ 2 C V(r) dr 

< C2e Jo ^ ^ , 

where C2 depends on t and n. Thus, for some C3 that depends on t and n, 

E[RlYt)]<C3(3nN-^E\ f 
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For the second remainder, we have the following estimate, using the 
Cauchy-Schwarz inequality: 


E[Rl2it)]<t^ [sup'finis) jE 


'■*^2 fjV{r)-\„)dr 


ds 


Uo 


< c{nYeH^E 




ds 


Hence, 

(3.11) 


E[{Rr,^Yt) + RndW] 

< 2(c3/3n^”^ +C(n)2e2t2)^[g2/Jy(r)drj^ 


Recall that, for some n, there exists a > 0 such that infAr^e |l^n(0)| = a. Let 
no be the smallest n satisfying this condition. It follows from (3.10) and 
(3.11) that 


(3.12) 


£;[«(t)-n„„(0)e/o(^M-^"o)'^’')2] 

< 2(c3/3no^^“^ + c(no)^e^t^)Ll[e^-^o 


Suppose that N is large enough (and, consequently, e is small) so that 

(3 13) ^ Ar—1 I ^2^2J.2^ ^ /I /o^„22A„r,t 

Since I 


4(c3/3no^ ^ + c(no)Vt^)<(l/2)a^e 


‘^no I loo 


< oo and e“ J2x l^(^) t)| = 2, we obtain 


(3.14) \unoit) \ — 

By (3.12), (3.13) and (3.14) 


e'^ u{x,t)4>noix) 

x&De 


< C4. 


^2g-2A„ot^[g2/>(r)drj 

<F;[(n„o(0)e-/'o('^W-^"o)'^")2] 

< 2E[«(t) - <(0)e/o (^W-^"o)'^-)2] + 2E[(n^Yi)f] 

< 4(c3/3nolV“^ + c(noYe^t^)E[e^fo + 2 C 4 

< (l/2)o2e"2^"oi£’[e2/o +2cl. 


so 


(3.15) 


^[g2/;y(r)dr-]<4^2^-2g2A„ot_ 
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Note that Ptu{0)^ = e~^'^^Un{0) and let v{x,t) = e-lo 
We combine (3.10), (3.11) and (3.15) to see that if —> oo (and, therefore, 

£ —f 0), then, for every fixed n and t > 0, Un{t) — Vn{t) —> 0 in distribution. 
Then Lemma 2.2 shows that 

(3.16) ^ eV=(a;,t)i,,-e/o^W'^^Ptu^-^(0)^5o in Mi(Mi.(;D)). 

x£De 

Step 3. It is not obvious that the normalization of Ptu^'^{0) in (3.16) is 
the same as in the statement of Theorem 3.2. It is conceivable that a sizeable 
proportion of positive and negative particles are tightly interspersed so that 
their masses cancel each other in the limit. We will show that this is not the 
case—intuitively, the two populations occupy disjoint parts of D. 

Let Bs{x) denote a hypercube in Hg, centered at x, with side length 6. 
We will consider only 6 > e. We set 

A5(Lx) = minl ^ % 

\y£Bs{x) y£Bs(x) 

Note that if e < 5/2 and x G D^, then Bs{x) contains at least ((5/2e)'^ sites 
y G Hg. Fix an arbitrarily small cq > 0. Suppose that, for some Bs{x), we 
have As{t,x) > co((5/2e)'^. Given this assumption, we will show that the 
number of (+ and —) particles that are located in B^{x) at time t and collide 
with a particle of the opposite sign before time t + <5^ has expectation greater 
than ciAs{t,x). Suppose, without loss of generality, that there are fewer + 
than — particles in Bs{x) at time t, that is, As{t,x) = J2yeBs{x) Vx ■ Consider 
independent continuous time reflecting random walks Yk, 1 < k < As{t,x), 
starting from the same points as the locations of the + particles at time t 
in Bg{x). The distribution of a process {yfc(s),s > t} is the same as for a 
single particle in our process rj, except that l^’s do not interact with other 
particles. For k< (l/2)((5/2e)'^ AAs{t,x), 

P{Yk{t + 5^) ^ Yj{t + 52), 1 < j < A; - 1) > Pi > 0, 

where pi depends only on the dimension d. This implies that, with probabil¬ 
ity p 2 , the number of different sites occupied by Yk{t + 5'^), 1 <k < As{t, x), 
is greater than C2(5/e)'^, where P 2 ,C 2 >0 depend only on d and cq. 

Choose As{t,x) locations of the — particles at time t in Bs{x). From each 
of these points, start a continuous time reflecting random walk Assume 
that Zfc’s, 1 < A; < As{t,x), are independent, and they are independent of 
Yfc’s. Suppose that the number of different sites occupied by Yi^{t + 5'^), 
1 < k < As{t,x), is greater than C2(5/e)'^ and call the set of these sites F. 
Then it is easy to see that, with probability greater than ps > 0, the number 
of distinct sites in T occupied by Zfc(A+ 5^), 1 < A; < (c2/2)(5/e)'^, is greater 
than 03 ( 5 /e)^. 
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Suppose that there are at least C2((5/e)'^ sites in F and the number of 
distinct sites in F occupied by Zk{t + 5^), 1 < A: < (c2/2)((5/e)'^, is greater 
than c^{5/eY. Find the first time ti> t when some and Zj occupy the 
same site and call these particles “eliminated.” Then, by induction, find 
the smallest tm > tm-i when some noneliminated Yk and Zj occupy the 
same site and eliminate this pair of particles. Note that the total number of 
eliminated pairs by the time f + 5^ cannot be smaller than {c^/2)[5/e)^. 

Now we return to our original model, with interactions between particles. 
Consider the set of + and — particles in rj that reside at the same locations 
as Tfc’s and Z^s at time t in B^{x). Choose from this set a pair ((5+,Q_) 
consisting of a + and a — particle and suppose that it would have been 
“eliminated” in the scheme described above, that is, if these two particles 
had been and Zj for some k and j. If (5+ and Q- do not meet before 
time f + (5^, it means that one of these particles must have met a particle of 
the opposite sign (different from Qj^ and Q-) before time f + and hence, 
at least one of particles <5+ and Q- has a jump before time t + 5^. Thus, 
with probability greater than P 2 P 3 > 0, there will be at least {c^/Y){5/eY 
jumps between times t and t + 6'^, by particles that are located in Bs{x) at 
time t, assuming that As{t,x) > co{6/2eY. 

Let K{t) be the number of collisions before t. Then EK{t) = {N/2)EJ^V{s)ds. 
Let Hs{t) = {x & Ds :As{t,x) > cq{ 5/2eY}. We see that, for some ci depend¬ 
ing on Co, 

EK{t)> E ^ E ci[As{k5\x)]. 

k=0 xeHs(kS^) 

For the same reason, for s € [0,(5^], 

EK{t) > E ^ E ci[Ai{k5^ + s,x)]. 

k=0 xGHs(kS^) 

rt-2S^ 

EK(t)> / 6~^E E ci[A5(s,a;)] ds. 

^ xGHs(s) 

For a function / on D^, let E{x,6,f) = J2yeBs(x)fiy)- Recall that ut = 
u{x,t) = We have 

(3.18) E{x,5, \ut\) - \E{x,6,ut)\ = 2N~^e~'^As{t,x). 

If X ^ E[s{s), then either 

E dy <co(d/2e)'^ or E dy < co(d/2e)'^- 

yGBsix) yGBsix) 


Hence, 

(3.17) 
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Hence, 


(3.19) 


(H(x,5, |ms|) - |H(x,5,Us)|) < Y co(5/2e)'^ 

x^Hs{s) x^Hs{s) 

< Y co{d/2s)^ <C 0 C 4 S~^, 

xGDs 


where C 4 is a constant depending only on D. Recall that N > c^e In view 
of (3.15), (3.17), (3.18) and (3.19), 


(3.20) 


rt-25'^ 

L 

= E 


2 - Y \-{x,S,£‘^Us)\ j ds 

^ x£Ds / 

f-t-2<52 


= £<^E 


= £‘^E 


/ Y {'^{x,5,\e‘^Us\)-\E{x,5,£'^Us)\)ds 
xeDs 

pt—25^ 

•152 


xGDs 


/ Y (S(a:,5, Im^I) - |H(x,(5,Ws)|)ds 


+ £^E 


x^Hsis) 
1-1-26'^ 

<52 

d 

l‘t-25'^ 

<52 

r-U2; 


< coC4t + 2e'^i? 


Y (H(x,5, |rts|) - |H(x,5,tt^)|)(is 

x&Hs(s) 


Y^ N ^£ ‘^As{s,x) ds 

xeHs(s) 


<coC4t + 2cQN ^5‘^EK{t) 

= CQCit + 2 cq 5^E [ V{s)ds 

Uo 

< CoC 4 t + C7(5^, 


where C 7 depends on cq. For a fixed t > 0 and n > 0, we can find cq > 0 and 
6 n G (0,2“") so small that the right-hand side of (3.20) is less than 2“". Let 
Tn be the set of s G [0,t] such that, for e < (5„, 


2- Y l^(x,d,£'^Us)l 

xeDg^ 


>n^2-^. 
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Then, by (3.20), \Tn\<n Choose arbitrarily small c* > 0 and let no be so 
large that Y.n>no < c*. Let T; = (0,t] \ [jn>no'^n- For and e < (5„, 

(3.21) pU- \^{x,6,e^Us)\>n^2-A <n-^. 

V / 


By passing to a subsequence, if necessary, we may assume that xeDe 
s)i 2 ; converges in Mi(Mf(D)). It follows from (3.21) that any limit of 
SxeDe i® supported on measures with the total variation 2, 

for every s €%,. Since c* > 0 in the definition of is arbitrarily small, we 
obtain the same conclusion for almost every s G [0,t]. This and (3.16) imply 
that, for almost every s > 0, 

(3.22) ^ s)i, - PsU^’^iO) ^ do in Mi{Mf(D)). 

xGDg 


Step 4. We will next show that, for large N, the process r —>■ V(s) ds is 

close to being continuous, in the sense that 


limsuplimsupPf sup 

<5i, 52—>0 e^O \0<ti,t2<t,\ti—t2\<Si 



V (s) ds > 62 


= 0 . 


We need the above claim in the last step of the proof. 

Let C{s) be such that C'(s)Pstt^’^(0) = PsU^’^{0). Note that C{s) G (0, 00 ) 
for all s because the solution to the heat equation in D with Neumann 
boundary conditions and nonzero initial condition with finite total variation 
is never identically zero and for all s, it has a finite variation. It is easy to 
see that (^(s) is a nondecreasing function. By assumption, (7(0) = 1. Recall 
that, for some n, there exists o > 0 such that infAr^e |lin(0)| = a and Psu{ 0 )^ = 

e~^^^Un{0). Since Psu{ 0 )^ = J 4>ndPsu{0) and (l)n is bounded (see the proof 
of Lemma 2.2), it follows that the total variation of Ptu{0) is bounded below 
and, therefore, C{t) is bounded above by a constant depending on D,n,a 
and t. 

By Theorem 2.1 of [3], for any 0 < si < S 2 < 00 , the Neumann heat kernel 
Ps{x, y) is Holder continuous jointly in (s, x, y) on [si, S 2 ] x Dx D. Let f{s, x) 
be the density of Psu{0) at x G H. We see that, for any 0 < si < S 2 < 00 , 
the family of densities {/(s,x)} corresponding to all measures sat¬ 

isfying our assumptions is equicontinuous on [si,S 2 ] x D. This implies that 
the family of functions {(7(s),0 < s <t} is equicontinuous. Since the fam¬ 
ily {/(s,x)} is equicontinuous and the total variation of Psu{ 0 ) is bounded 
below on [si,S 2 ] x D, it follows that the Prohorov distance of Psu{0) from 
the measure identically equal to 0 is bounded below on [si, S 2 ] by a constant 
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depending only on D,n,a,si and S 2 - By (3.16) and (3.22), for almost every 
s > 0, 

(C(s) - e/o ^ (5o in Mi{Mf(D)). 

We see that C{s) — e-^o must converge to 0 in distribution for al¬ 

most every s. Since C{s) are uniformly continuous and nondecreasing, and 

e/o is nondecreasing, it is easy to see that C{s) - e/o converges 

to 0 uniformly on compact intervals, in distribution, as —> oo. This implies 
the claim stated at the beginning of Step 4. 

Step 5. Next we will show that the convergence holds not only for every 
fixed t > 0, but also in the Skorohod topology on Mi{S{{0,oo), Mf{D))). 
Fix a smooth function ip on D and let = wt = {ut, ip). We will show that 
the family of processes {w^.,£ > 0} is tight in Mi(5((0,oo),R)). In order to 
prove that, we will first derive some estimates for wt similar to the estimates 
for Un- We have 

(3.23) dwt = Lwtdt + dM^(t), 

where M^(t) is a martingale. Since ip is smooth, there exists a series expan¬ 
sion for ip that yields 

^eV{x) = A(/?-L V'e(x), 

where |V'e(3^)| < We have Lrjx = A*rya; -|- Vrjx and 

Wt = {u,ip) =£'^ N~^£~‘^r]xip{x) = N~^ ^ Vx^{x), 

x£De x£Ds 

SO 

Lw = N~^ ^ {Lr]x)ip{x) 

x&De 


= N ^ {/Slr]x)ip{x) + N ^ Y i^dx)^{x) 


x&De 


x£Ds 


= ^ ^ ^iy)Pyi>idy - K ^{x)PxyVx) I (fix) 

x&De \yeDe ) 


KyeD, 

+ N~^V Y dxip{x 


xeDe 


= ^ ^ Z ( Z iK^ix)pxyip{y) - h^^{x)pxyip{x))\r]x 

xeDe \yeDp_ } 


-LA V ^ r]xip{x) 

x&De 


= A ^ Z rjxi^e^pix)) + VN ^ Y Vx<f{x) 

x£Ds x£De 


(3.24) 
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= N ^ r]a;{A(p{x)+ipi;{x))+ VN ^ Y Vx^ix) 

xGDs x^De 

= {u,A(p) + ^^ + Vw, 
where l^'el < c^e. Note that 

(3.25) E[M^it)]=J^E[A^{s)]ds, 

where = A^p{t,rj) is given by the following formula: 

K = lim(l/s)£;[(w;(t + s) - w{t)f\r]{t)] 

s—>0 

= iV-2 Y K^{x)Pxy\'n+l{^^yQ^{^{y)-^{x)f 

x,veDp_ I 


+ vt^{riy<0}N~‘^ Y vtVv{piy)-pix) 

u^vGDs 

+ (f{u) - (p{v)y 

+vy^{yy>o}^~‘^ Y yyvtiv’ix)-(p{y) 

u,vGDe 

- (f{u) + (p{v)y 


We obtain the following bound for A^^ for small e: 

(3.26) < c,N-^V{ey\V^\\l + ||:/.||L) < ^N-^V, 

where (5 <oo depends on tp. 

It follows from (3.23) and (3.24) that 


e ^0 


(3.27) 


w{t) = eIo^^^^'^^(w{0) + 

+ J e~-l^o^^'"'>^'~({u,Ap){s) + ^eis))ds^ 

w(0)efo efs dM^(s) 

Jo 

+ [ ds + [ {u, Ap){s) ds 

Jo Jo 


df 

= W 


(O)efo (t) + R 2 (t) + Ro(t). 
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In view of (3.25), (3.26) and (3.15), we have 
E[Riit)]=E 




< C2PN-^E 


Uo 


= C2/3iV^£;[(l/2)(e^-/'o - 1)] 

<C3/3N-^Ele^fo^(’'^'^^] 

< C4e^^"o‘iV“^ 

Since i?i(f) = Jq dM^{s) and M^{t) is a martingale, so is i?i(t). 

Hence, by Doob’s inequality. 


E 


sup Rl{s) 
L0<s<t 


< 4£;[ii^(t)] < 4c4e2^"o*iV“^ 


A similar calculation and (3.15) show that 


E 

Hence, 

(3.28) 


sup R^is) 

0<s<t 


f sup '^‘^{s)]e [ ds 

\o<s<t ) Uo 




E 


sup (i?l(s) + R 2 {s)y 

0<s<t 


<cg{N ^ + eV), 


where cg depends on t and ip. 

For the last term on the right-hand side of (3.27), we observe that, for 
0 < fi < f2 < 


|/?3(f2) - R3{ti)\ = 


C^2 rt‘ 


'h 


y A(y9)(s) ds 


+ f ‘^^){u,Aip)is) ds 

Jo 

f e-^s {u, Aip){s) ds 


< 


+ [ ^ els' V(r),r 

Jo 

sup e-^s f |(tt, A(y9)(s)| ds 

i<s<t Jtl 

I f IN fW(r)dr ,/ A N| , 

+ (e-'*! —1) sup e-ix ^ ' / |(tt, A(y9)(s)| ds 

0<s<t Jo 
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OO • 


In view of (3.15) and the first claim in Step 4, 

lim sup lim sup P ( sup | Rs {t 2 ) — Rsiti 

<51,52^0 e^O \o<ti,t2<t,\ti—t2\<Si 


■0<ti ,t2<t,|ii—t2j<<5i 


sup \R3{t2) - R 3 {tl)\> 62 ) =0. 


This, (3.27), (3.28) and Lemma 2.3 show that {w‘^,e > 0} is a tight fam¬ 
ily of processes in Mi(5((0,oo),M)). Since smooth functions are dense in 
the set of continuous functions on D and the sum of two smooth func¬ 
tions is smooth. Theorem 3.7.1 of [12] shows that the family of processes 
{SxgDe Oix, e > 0} is tight in Mi(5((0, oo), Mi7’(L)))). Since 

{Po'^^'^’'^(0), e > 0} are tight and the process t Ptu^’^{0) is a continuous 
function with values in Mp{D), completely determined by Pqu^'^{Qi), we 
conclude that 



is tight in Mi(S'((0, oo), Mi?(Zl))). By (2.13), any convergent subsequence of 
this family must be identically equal to Joj by the right continuity. □ 

Acknowledgments. We are grateful to Rodrigo Banuelos, Victor Ivrii, 
John Sylvester, Alain Sznitman, John Toth and Steve Zelditch for very help¬ 
ful advice. We thank the referees for many helpful suggestions. 


REFERENCES 


[1] Banuelos, R. and Burdzy, K. (1999). On the “hot spots” conjecture of J. Rauch. 

J. Func. Anal. 164 1-33. MR1694534 

[2] Banuelos, R. and Pang, M. (2006). Level sets of Neumann eigenfunctions. Indiana 

Univ. Math. J. 55 923-939. MR2244591 

[3] Burdzy, K. and Chen, Z. (2002). Coalescence of synchronous couplings. Probab. 

Theory Related Fields 123 553-578. MR1921013 

[4] Burdzy, K., Holyst, R., Ingerman, D. and Margh, P. (1996). Configurational 

transition in a Fleming-Viot-type model and probabilistic interpretation of 
Laplacian eigenfunctions. J. Phys. A 29 2633-2642. 

[5] Burdzy, K., Holyst, R. and March, P. (2000). A Fleming-Viot particle repre¬ 

sentation of Dirichlet Laplacian. Comm. Math. Phys. 214 679-703. MR1800866 

[6] Conti, M., Terracini, S. and Verzini, G. (2005). Asymptotic estimates for the 

spatial segregation of competitive systems. Adv. Math. 195 524-560. MR2146353 

[7] Conti, M., Terracini, S. and Verzini, G. (2005). A variational problem for the 

spatial segregation of reaction-diffusion systems. Indiana Univ. Math. J. 54 779- 
815. MR2151234 

[8] Courant, R. and Hilbert, D. (1953). Methods of Mathematical Physics. Inter¬ 

science Publishers, New York. MR0065391 




AN ANNIHILATING-BRANCHING PARTIGLE MODEL 


25 


[9] Cybulski, O., Babin, V. and Holyst, R. (2004). Minimization of the Renyi entropy 
production in the stationary states of the Brownian process with matched death 
and birth rates. Phys. Rev. E 69 016110. 

[10] Cybulski, O., Babin, V. and Holyst, R. (2005). Minimization of the Renyi en¬ 

tropy production in the space partitioning process. Phys. Rev. E 71 046130. 
MR2139992 

[11] Cybulski, O., Matysiak, D., Babin, V. and Holyst, R. (2005). Pattern formation 

in nonextensive thermodynamics: Selection criterion based on the Renyi entropy 
production. J. Chem. Phys. 122 174105. 

[12] Dawson, D. A. (1993). Measure-valued Markov processes. Ecole d’Ete de Probabilites 

de Saint-Flour XXI — 1991 1-260. Lecture Notes in Math. 1541. Springer, Berlin. 
MR1242575 

[13] Ethier, S. N. and Kurtz, T. G. (1986). Markov Processes. Characterization and 

Convergence. Wiley, New York. MR0838085 

[14] Grieser, D. (2002). Uniform bounds for eigenfunctions of the Laplacian on manifolds 

with boundary Comm. Partial Diff. Eg. 27 1283-1299. MR1924468 

[15] Grigorescu, I. and Kang, M. (2004). Hydrodynamic limit for a Fleming-Viot type 

system. Stochastic Process. Appl. 110 111-143. MR2052139 

[16] Grigoresgu, I. and Kang, M. (2006). Tagged particle limit for a Fleming-Viot type 

system. Electron. J. Probab. 11 311-331. MR2217819 

[17] Hempel, R., Seco, L. A. and Simon, B. (1991). The essential spectrum of Neumann 

Laplacians on some bounded singular domains. J. Func. Anal. 102 448-483. 
MR1140635 

[18] Morrey, C. (1966). Multiple Integrals in the Calculus of Variations. Springer, New 

York. MR0202511 

[19] Netrusov, Yu. and Safarov, Yu. (2005). Weyl asymptotic formula for the Lapla¬ 

cian on domains with rough boundaries. Comm. Math. Phys. 253 481-509. 
MR2140257 

[20] Protter, Ph. (1990). Stochastic Integration and Differential Equations. A New Ap¬ 

proach. Springer, Berlin. MR1037262 

[21] Stroock, D. and Varadhan, S. R. S. (1971). Diffusion processes with boundary 

conditions. Comm. Pure Appl. Math. 24 147-225. MR0277037 


Department of Mathematics 
University of Washington 
B ox 354350 

Seattle, Washington 98115-4350 
USA 

E-mail: burdzy@math.washington.edu 


Department of Mathematics 
University oe Toronto 
100 St. George St. 

Toronto. Ontario 

Canada MSS 3G3 

E-mail: quastel@math.toronto.edu 


